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ABSTRACT 

From Sloan Digital Sky Survey u'g'r'i'z' imaging, we estimate the stellar masses of the host galaxies 
of 70 low redshift SN la (0.015 < z < 0.08) from the hosts' absolute luminosities and mass-to-light 
ratios. These nearby SN were discovered largely by searches targeting luminous galaxies, and we 
find that their host galaxies are substantially more massive than the hosts of SN discovered by the 
flux-limited Supernova Legacy Survey. Testing four separate light curve fitters, we detect ~2.5er 
correlations of Hubble residuals with both host galaxy size and stellar mass, such that SN la occurring 
in physically larger, more massive hosts are ~10% brighter after light curve correction. The Hubble 
residual is the deviation of the inferred distance modulus to the SN, calculated from its apparent 
luminosity and light curve properties, away from the expected value at the SN redshift. Marginalizing 
over linear trends in Hubble residuals with light curve parameters shows that the correlations cannot 
be attributed to a light curve-dependent calibration error. Combining 180 higher-redshift ESSENCE, 
SNLS, and HigherZ SN with 30 nearby SN whose host masses are less than 10 10 8 M in a cosmology 



fit yields 1 + w = 0.22_ 01 q 8 , while a combination where the 30 nearby SN instead have host masses 
greater than 10 10 8 M Q yields 1 + w = -0.03^;^. Progenitor metallicity, stellar population age, and 
dust extinction correlate with galaxy mass and may be responsible for these systematic effects. Host 
galaxy measurements will yield improved distances to SN la. 
Subject headings: supernovae: general 



1. INTRODUCTION 



(Astier et al.||2006; Wood-Vasey et al.||2007| |Riess et al. 


2007 Hicken et ai.| 2009 b 


Kessler et al.||20U9a). Their 



tive 56 Ni during the explosion and subsequent decays 
to 56 C o and 56 Fe ( [Truran et al.||l967| [Colgate fe McKee 
1969 1 then power the SIN light curve. How the progenitor 



shifts and, after calibration by light curve shape, yield 
luminosity distances over a wide range of cosmic epochs 
spanning the past 8 Gyr. While the peak luminosities 
of SN la vary by a factor of ^3, more slow ly-evolving 
SN la are intrinsi cally brighter ([Phillips 1993) and bluer 
( Riess et al.|1996 ), a pattern used to calibrate SN la and 
measure luminosity distances to a precision of ~10%. Al- 
though initial efforts were limit ed by small detector sizes 



accumulates a mass close to the Chandrasekhar limit be 
fore the explosion is not yet clear, but plausib le scenarios 
include accretion from a binary comp anion (Whelan & 
Iben|1973| |Han fc Podsiadlowski|200|| _ o r a m e rger with 
a sec ond white dwarf (Iben fc Tutukctv||l984| |Webbink 



( Norgaard-Nielsen et al. 1989), search programs using 



newer, large format cameras discovered that distant SN 
la were ~25% fainter than expected for any matter-only 
universe, favoring an accelera ting expansion ( |Riess et al. 
1998[ |Perlmutter et al.|l999 1 . In this paper, we look for 
correlations between the Bubble residuals of nearby SN 
la and their host galaxy sizes and masses. 

Theoretical models of a carbon-oxygen white dwarf 
with a mass approaching the Chandrasekhar limit show 
that it can explode to produce a rea li stic SN la spec- 
trum jHillebrandt fc Niemeyer| |2000l |Kasen fc Plewa 
[2005] [Kasen fc Woosley][2007| |Kasen et al.j|20o9| ). A su- 
personic deflagration wave, possibly preceded by a more 
slowly moving carbon fusion flame, synthesizes radioac- 



1984J). 

Galaxies with little or no star formation can host SN 
la, implyi ng that some fraction of SN la have old pro- 
genitors. IQemler fc Tinsley|[T979"l showed that the SN 
la rate in star-forming galaxies is related to the current 
rate of star formation, pointing to a second, rapid path to 
SN la. Subsequent studies have confirmed an increased 
SN la rate per u nit mass in galaxies with younger stel- 



lar p opulations (van den Bergh 1990 Mannucci et al 



20051, which has been used to rule out strongly a sim- 



ple progenitor population co nsisting of stars that explode 
after reach ing a single age ([ Sullivan et al.||2006| |Prieto 
et al.|200"8j ). Progenitor populations with detonation de- 
lay times in cluding two separate compon ents, 'prompt' 



and 'tardy' ( Scannapieco fc Bildsten|2005 ), or that form 
a continuum better fit the SIN la rate. Consistent with 
a relatively young progenitor population in star-forming 
hosts, the distribution of SN la in the g '-band light of 
spiral galaxies is similar to that of S N II (Kolmogorov 
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Smirnov p=0.22; |Kelly et al.||2008[ ) 
find that core collapse SIN (including SJN Ib/c and S" 
are more closely associated with brighter (/'-band pixels 
than SN la when the regions outside of the bulge are 
considered in grand-design spirals (Kolmogorov-Smirnov 
p=0.004), although whether SN Ib/c should be included 
in this comparison is not clear given their preference for 
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chemically enriched environments ( Prieto et al.|2008 ) of- 
ten found in inner, higher surface- brightness regions. 

Correlations between Hubble residuals and progenitor 
properties have the potential to introduce bias into SN 
la cosmological measurements, because SN at high red- 
shift are likely to have younger and possibly less metal 
rich progenitors than their low redshift counterparts. If 
correlations are identified, host measurements such as 
size, mass, or metallicity could be combined with light 
curve para meters in, for example , a hierarchical Bayesian 
framework ( Mandel et al.| [2009) to improve luminosity 
distance estimates to SiN la and control systematics in 
cosmological measurements. As of yet, however, there is 
no consensus on whether any correlations between host 
galaxy properties and SN la Hubbl e residuals exist 



(Jha et al. 2007) with Ry 

(MLusaiJ. — 



1.7 (MLCS17) and Ry = 3.1 



From spectra of 29 E/SO galaxies, | Gallagher et al. 2008 
found that host metallicity correlates with MLCS2k2 
Hubble residuals at the 98% confidence level for a one- 
sided test. For a two-sided test, in which either po sitive 
or negative correla tions would be interesting, the Gal- 
lagher et al. 2008 confidence would instead be ~96%, 



a 'two-sided' confidence. Unlike Galla gher et al.||2008| 
IHowell e"t~a l. 2009 did not find a trend in Hubble residu- 
als of 57 SN from the flux-limited SNLS lst-year survey 
using SiFTO to fit light curves. These studies did not at- 
tempt to separate host dependence from light curve de- 
pendent trends in Hubble residuals, important given the 
correlatio n between light cu rve shape and host galaxy 
type (i.e. |Hamuy et al . |[l996] ) . 

Here we look lor correlations between the Hubble resid- 
uals of nearby (0.015 < z < 0.08) SN la and host galaxy 
properties using four separate light curve fitters. We si- 
multaneously fit for trends in Hubble residuals with both 
light curve and host galaxy properties to isolate a de- 
pendence of the SN la width-color-luminosity relation 
on host environment. In Section [2] we briefly discuss 
the light curve fitters used to analyze SN in this paper. 
Section [3] reviews current understanding of the effect of 
progenitor properties on SN la luminosities and Hubble 
residuals. The data set and methods used to measure 
host galaxies sizes, masses and colors are described in 
Section|4] We present and examine the significance of the 
trends we find in Hubble residuals with light curve pa- 
rameters, host galaxy size, and host galaxy stellar masses 
in Section [6] Finally, in Section [7| we discuss and inter- 
pret these trends. In this paper, H o =70 km s _1 Mpc -1 is 
used to calculate galaxy sizes and st ellar masses. A con 



cordance cosmology is assumed (see |Hicken et al. 2009b 



for parameter values inferred from cosmology tits) 
2. LIGHT CURVE ANALYSIS 

We use the light curve fits and Hubble residuals from 
the |Hicken et al.|[2009b| cosmology analysis. Their low 
redshift SN sample was constructed from ~70 existing 
light curves in the literature and ~13 new light curves 
(the CfA3 sample; [Hicken et al.|2009a| ). CfA3 light curve 
observations were acquired with the 1.2m telescope at 
the F. L. Whipple Observatory and reduced using a sin- 
gle pipeli ne. A high percentage of these SN have spec- 
troscopic (Matheson et al.||2008| Blondin et al. 2010, in 



preparation) and near-infrared photometric data ( Wood 
Vasey et aT] 2008[ Friedman et al. 2010, in preparat ion 
Fitters app lied to the light curves inc lude: SALT (Guy 
etaT1[2505| ), SALT2 ( [Guy et ai.||2"007) , and MLCS2k7 



2.1. MLCS17 and ML CS31 



In the MLCS2k2 ( |Jha et al.||2007[ ) model, the light 
curve color, the light curve shape, and the SN intrinsic 
luminosity are a function of one parameter, A. SN that 
have light curves with low A decay more slowly and are 
intrinsically brighter and bluer. MLCS2k2 assumes that 
all SN la follow this A parameterization of light curve 
width, color, and SN luminosity and that any remaining 
variation arises from extinction by dust, which both di- 
minishes the apparent brightness of a SN and reddens its 
color according to the relationship Ry = Ay /E(B — V) 
where Ay is the extinction and E{B — V) is the color 
excess. Ry parameterizes the wavelength dependence of 
the dust reddening law. To fit a light curve, spectral SN 
la templates are first transformed to the observer-frame 
(as opposed to SN r est-frame) through K -corrections us- 
ing the procedure in |Nugent et al. ||2002| MLCS2k2 then 
simultaneously fits for A and the extinction Ay using a 
one-sided exponential prior with scale te(b-v) = 0.138 
for low reds hift SN. For higher redshift SN, |Hicken| 
|et al.|[2009b| uses the 'glosz' prior, which incorporates 
a detection effi ciency that depends on redshift, for the 
ESSENCE SN (|Wood- Vasey et al.|[2fWL and a modified 
'glosz' for the Higher-z SIN ( |Riess et al.||2005J ). 

While the mean Milky Way Ry value is 3.1, this value 
may not accurately describe extinction in other galaxies. 
|Astier et al |2006| and |Conley et al 12007] find that 0, the 
coefficient between luminosity and SN color, is close to 
2, less than the 4.1 expected if Ry = 3.1 and the SN 
la color variation is due entirely to dust. By compar- 
ing average colors of SN la at a series of epochs against 
simulated co lors inside the MLCS2k2 framework, |Kessler] 
|et al.||2009a| measured Ry from the SDSS-II SN sample, 
finding that Ry = 2.18 ±0.14(stat)±q.48(sys t). Recog- 
nizing this uncertainty, Hicken et al. 2009b performed 
fitting using both Ry = 3.1 and Ry = 1.7, the value 
that minimizes the scatter in the Hubble diagram, la- 
beling these as MLCS31 and MLCS17. The MLCS2k2 
model parameters were fit from many high-quality, low 
redshift SN including a number of 1991bg-like SN, which 
are a faint class of S N la identified by strong Ti II lines. 



Hicken et al. 2009b note that removing the 1991bg-likc 
SN from MLCS2k2 training sets may yield i mproved cal- 



ibratio n of 'normal' SN la [see Figure 26 of |Hicken et al. 
2009al. 



2.2. SALT and SALT2 

SALT and SALT2 fit for a stretch parameter (s for 
SALT, x x for SALT2), a SN color, c = [B - V) t=Bm ax + 
0.057, and the time of maximum _B-band light, t$. High 
stretch SN are more intrinsically luminous and have 
broad, slowly-decaying light curves. SALT uses a mod- 
ified Nugent template ( 2002 1 to fit each SN light curve 



and derive light curve parameters. A linear function of 
s and c gives the distance modulus in the SALT model, 
including the effects of intrinsic color variation and dust 
reddening in the single c term, 



fl B — m B 



M + a(s - 1) - /3c 



(1) 



where the coefficients a, ft, and M are marginalized over 
in the cosmological fit. Instead of using the Nugent tern- 
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plate, SALT2 builds a light curve model as a function of 
stretch and color from a combination of empirical light 
curves and spectra from low and high redshifts. This 
approach benefits from observations of high redshift SN 
that sample the rest frame UV and creates a model rep- 
resentative of SN light curves at low and high redshift. 
The stretch term is x\ and the color c, 

[iB = m^ ax — M + axi - /3c (2) 

where the linear coefficients a, (3, and M are also 
margin alized over in the cosmological fit. |Hicken et al.| 
|2009b| found a = 1.26 and /3 = 2.87 in the case of SALT 
and a = 0.109 and p = 2.67 for SALT2. SALT and 
SALT2 do not include 1991-bg like objects in training. 

3. PROGENITOR PROPERTIES AND SN IA LUMINOSITIES 
3.1. Trends in Total SN Luminosities 



may alter the SN la width-color-luminosity relati on suf- 



|Hamuy et al. 1996 noticed that SN la with fast- 
declining light curves and low luminosities are more com- 
mon in early-type than late-type hosts, the first strong 
indication that progenitor properties influence the explo- 
sion energy. Relationships have been found between SN 
la luminosity and the star form ation rate in a flux- li mited 
survey from host photometry (Sullivan et al.|2006|), stel- 



lar population age from sp ectroscopy ([Gallagher et al 



2005 [Gallagher et al. 2008[), and stellar pop ulation age 
from host photometry ( [Howell et al. 2009). However, 



these studies find similarly strong trends in SN la lumi- 
nosity with host metallicity, leaving open the questi on of 
whether stellar age or metallicity is mo re importa nt (|Gal- | 
lagher et al.l|2008| |Howell et al.||2009[ ). Recently, |HicEen| 



et al.||2U09b| found that the brightest SN la occur m in- 
termediate spirals, not the galaxies with the youngest 
stellar populations. 

Galaxies have metallicity gradients such that stars far- 
ther from the galaxy center are generally less metal-rich, 
and studies have exami ned SN la properties at different 



projected radial o ffsets (Wang ct al. 1997 [ Ivano 1 
Gall agher et al. 2005| Boissier &; Prantzos 
1997[and |Glallagher et al.|iW05| found! 



Ivano v et ah 
2009[) . 
fftaTSN 



2000 

Wang et al.| 

la at larger galactocentric radii have less luminous explo 
sions. However, it is likely that this trend would disap- 
pear if SN la hosted by large elliptical galaxies, which 
account for many of the SN la at large galactocentric 
radii, were removed. Large ellipticals host less luminous 
SN la at both sma ll and large galactocentric distances. 
|Cooper et al.|2 009 recently studied the local galaxy den- 
sity, the number of galaxies per unit volume, near the 
hosts of SDSS-II SN la. They found that, when SN la 
were detected in blue, star-forming galaxies, a dispropor- 
tionate percentage of the host galaxies were in regions of 
low galaxy density, possibly suggesting that prompt SN 
la are more luminous or occur more frequently in envi- 
ronments with lower gas-phase metallicities. They found 
no evidence for an analogous effect among SN la with 
redder, older host galaxies. 

3.2. Trends in Hubble Residuals 

|Timmes et al. |2003| suggested a plausible mechanism 
by which prog enitor metallicity c ould affect SN luminos- 
ity. Although Howell et al. 2009 found that variation in 
SN la luminosities predicted by the theory could only 
account for 7-10% of the observed variation, the effect 



ficiently to p roduce trends in Hub ble residuals ( Kasen 



et al. 20091. Timmes et al. 2003 argue that metal in- 



creases the nucleon density in the white dwarf before the 
explosion, and greater nucleon densities favor the synthe- 
sis of neutron-rich, iron-peak elements such as 58 Ni and 
54 Fe at the expense of 56 Ni during the burning phase. 
The diminished 56 Ni mass that results produces a less 

lumin ous explosion. 

The |Gallagher et"aLl[2005l and |Gallagher et~aLl [20081 
spectroscopic analyses of host galaxy properties con- 
strained trends in Hubble residuals with host galaxy 
metallicity using MLCS2k2 light curve fits. | Gallagher] 
et al. 1 120051 analyzed the integrated host spectra of 57 



SN la taken by panning a spectrograph slit across each 
galaxy. For the 28 galaxies in their sample with on- 
going star formation, they calculated metallicities from 
emission-line diagnostics, finding no strong evidence for 
a trend in Hubble r esiduals with gas-pha se metallicity. 
In a separate study, |Gallagher et al. 2008 measured the 
metallicities and ages of the stellar populations in E/S0 
hosts from absorption features, finding a correlation be- 
tween metallicity and Hubble residuals with a one-sided 
c onfidence level of 9 8 % . 

IHowell et al.l2009l constrained host metallicities less di- 
rectly by exploiting the correlation between galaxy mass 
and m etallicity (e.g. Tremonti et al. 12004 Gallazzi et al 



[2005J. They fit ugriz photometry with the PEGASE2 
population synthesis code to find host masses for SN 
from the flux-limited Supernova Legacy Survey (SNLS) 
and extrapolate d metallicities from st ellar mass measure- 
ments using the |Tremonti et al~|2 004 relation. They find 
no statistically significant trend in Hubble residuals with 
host metallicity, but th eir result may be consistent with 
|Gallagher et al. 2008| (the original paper reversed the 



reported direction ot trend in Hubble r esiduals; J. Gal- 



laghe r, private communication). Unl ike |Gallagher et al 



2008] who used the MLCS2k2 fitter, |Howell et al.pOW) 



ana lyzed light curves w ith SiFTO, a fitter implemented 
by ( |Conley et al.||20 08), who obtained relatively similar 
cosmological best-fit parameters from SiFTO and SALT2 
a nalyses of SNLS da t a . 
Neill et al.||2009] fit photometry of 168 low redshift 



SN host galaxies with the PEGASE2 population syn- 



thesis models used by |Howell et al.|[2009 For each host 
galaxy, they compiled GALEX UV and either Sloan Dig- 
ital Sky Survey (SDSS) u'g'r'i'z' or RC3 Johnson UBV 
magnitudes. In an effort to select SN la with low extinc- 
tion, they identified host galaxies that were best-fit by 
a PEGASE2 model with E[B - V) < 0.05. Of the SN 
hosted by these galaxies, they selected 22 with SiFTO 
light curve parameters s < 0.75 and c < 0.7 to produce a 
sample with light curve parameters similar to those of in- 
termediate redshift SNLS SN. Among these 22 SN, they 
found a 2.1a trend in SiFTO Hubble residuals with host 
galaxy age. However, they found no correlation between 
the PEGASE2 extinction and SN la color, so the trend 
is difficult to interpret. When all SN with s > 0.75 and 
c > 0.7 were fit together, no correlation between Hubble 
residuals and host age was seen. 

Galaxy morphology is another host property that may 
correlate with Hubble residuals. After correction for 



light curve shape and reddening, Jha et al. 2007 found 
that SN la in elliptical hosts were la brighter after light 
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curve correct ion with MLCS2k2 distances, while [Hicken| 
et al.||2009b| similarly found that SN la in E/SO hosts 
are brighter after light curve correction than those in 
Scd/Sd/Irr hosts by 2a. 



4. DATA 



The SPSS (York et al 



2000) DR 7 (Abazajian & 



Sloan Digital Sky Survey 2009) includes 11663 square 
degrees ot u'g'r'i'z' imaging taken in 53.9 second expo- 
sures with the 2.5 meter telescope in Apache Point, New 
Mexico. The SDSS filters cover the near-infrared detec- 
tor sensitivity limit to the ultraviolet atmospheric cutoff 
( Fukugita et al.|[l996 ). Each Sloan frame is a 13.5 x 9.9 
arcminute held over a 2048 x 1498 array of 0.396" pixels. 

4.1. Sample Selection 
Table [l] describes the selecti on of the sample from SN 



whose light curve fits meet the[Hickcn et al. 2009b 



imal cut' criteria for at least one light curve fitter. The 
'minimal cuts' are t lst < +10 days and xt < 1-6 for 
MLCS17, and t ut < +10 days as well as xl < U> for 
MLCS31. For the SAL T sample, |Hicken et al.||2009b| ap- 
plied the cuts used by |Kowalski et al.||2008| to construct 
the Union compilation including z > 0.015, t\ s t/s < +6 
days, fit convergence, a Hubble residual within 3er of 
best-fit cosmology, and five or more separate observa- 
tions. The SALT2 cut consists of excluding SN with 
light curves where xt > 10. 

From these SN, we select nearby (0.015 < z < 0.08) 
SN la in the Hubble flow whose host galaxy images were 
in the SDSS DR7. A further cut eliminates SN for which 
all SDSS host galaxy images w ere contaminated by SN 
light. In late-time light curves (Lair et al. 2006), SN la 
decay to fainter than ~-12 Johnson V magnitudes within 
one year. To avoid contamination from rising or declin- 
ing SN, we exclude images of SN hosts taken from three 
months before to one year after SN detection. We finally 
select SN without a nearby bright star and those with 
a MLCS17 light curve fit consistent with low extinction 
A v < 0.5 (R v = 1.7). 

For the MLCS17 and MLCS31 subsamples, we further 
restrict SN to those with A < 0.7, becaus e 0.7 < A < 
1.2 SN were found by Hickcn ct al. 2009b] to be poorly 
calibrated. 60 MLCS17, 58 MLCS31, 53 SALT, and 62 
SALT2 SN remain after all of the cuts. 

5. METHODS 

5.1. Photometry 

To determine SN la host galaxy colors and stellar 
masses, we measure and model host u'g'r'i'z' magni- 
tudes. SDSS photometry of bright galaxies suffers from 
overestimated local backgrounds and occasional galaxy 
'shre dding' by the object detection algorithm [see Fig. 



9 of Blanton et al. 2005 , so we perform new photom- 



etry measurements on SDSS images. To measure the 
total galaxy magnitude, we construct a mosaic of nearby 
frames. A mosaic includes the entire extent of the host 
light distribution so that the host magnitude can be mea- 
sured by directly summing photons. To construct each 
mosaic, we resample SDSS frames to a common grid us- 
ing SWarrj^Jwith the LANCZOS3 kernel, rescaling fluxes 

4 http://astromatic.iap.fr/ 



TABLE 1 

Sample Selection 



Criterion 


la 


Hickcn ct al. 2009b 


309 


U.Ulb < z < U.U8 


170 


SDSS DR 7 


89 


Residual SN Light 


79 


Bright Stellar Contamination 


76 


Low Extinction A v < 0.5 (Rv = 1.7) 


70 



Note. — Number of SN la remaining after 
apply ing each criterion. (1 ) SN with light curves 
fit by |Hicken et all] [2009b passing the 'minimal 
cuts'; (2) hosts with 0.015 < z < 0.08; (3) in- 
side SDSS DR 7 coverage; (4) hosts observed 
when expected residual SN light is insignificant; 
(5) no bright stars superimposed on galaxy or 
very bright stars nearby (removes SN 2003W, 
SN 2006bw, and SN 2007O); (6) objects with 
low extinction. 

according to SDSS zeropoints for each frame. Next, we 
determine background differences between the adjacent 
chips by measuring their median difference in their region 
of overlap. Subtracting away the median, we produce a 
constant background a cross the mosaic. 



With SExtractor's ( |Bertin fc Arnouts| [1996| dual- 
detection mode, we measure the photometric aperture 
MAG_AUTO on the g '-band image and apply this aper- 
ture to the u'g'r'i'z' pixel-matched mosaic images to 
measure total galaxy magnitudes on the SDSS data. Be- 
cause a consistent aperture is applied, galaxy light distri- 
butions in lower S/N images will not be inappropriately 
'shredded.' Low values of the contrast parameter DE- 
TECTJV1INCONT correspond to aggressive object de- 
blending. We use DETECTJV1INCONT=0.005 to ex- 
tract host galaxies except for the host of SN 2001en, SN 
2002jy, and SN 2007bz, where we use 0.03, and the edge- 
on host of SN 2001V where we use 0.5. MAG_AUTO 



uses a 



Kron ( 1980 ) radius, and our 'total' magnitudes 



sum flux within 2.5 Kron radii. The background level 
of each SDSS mosaic is estimated from the median of 
polygonal regions placed by hand outside the periphery 
of each galaxy Gala ctic extinction values are taken from 



Schlegel et al 



I Blanton et al 



1998, and we use KCORRECT (v. 4.13) 
2003 ) to calculate rest-frame magnitudes. 



5.2. Galaxy Masses 
5.2.1. PEGASE2 



With the code LePhare (Arnouts et al. 1999 Ilbert 



et al.||2006[ ), we fit our u'g'r'i'z' host magnitudes with 
spectral energy distributions (SEDs) f rom PEGASE2 
(Fioc & Rocca-Volmerange 1997 1999) stellar popula- 



tion synthesis models. The SEDs include both nebular 
and stellar light and track metal enrichment and metal- 
dependent extinction through 69 time-steps for each of 
the size mode ls. We use the initial mass function of lRanal 
|fc Basu||1992| with a 5% fraction of close binaries and a 
mass range of 0.09 to 120 M . Table [2] lists the star 
formation histories used to construct the SEDs. From 
the marginal probability distribution calculated from the 
likelihoods of SED fits to galaxy photometry, we find the 
median stellar mass and uncertainty around the median. 

Near-IR fluxes, which are not included in our fits, 
have greater sensitivity than optical measurements to 
the older stellar populations that account for much of 
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mean of residuals 
rms of residuals 



0.033 
0.15 



TABLE 2 

PEGASE2 Star Formation History Scenarios. 



v infall gal. winds 



extinction 



o 




-1 -0.5 

log(M PH0T PECASE2 /M D4000 Hs/ 

I I I I I I L_ 

11 



log M PH0TPEGASE2 
(a) 

Fig. 1. — Comparison bet ween our stellar masses fr om PEGASE2 
fits to u'g'r'i'z' fluxes and Kauffmann ct al. 2003 stellar masses 
that use spectral indices stellar Balmer absorption to constrain star 
formation history. Despit e using only u'g'r'i'z' SPSS photometry 
(petroMag) , we reproduce Kauff mann et aL| 2003 masses accurately 
with root-mean-square residuals ot U.lt> dex arid a mean bias of 
0.033 dex. Plotted are 10000 r andomly selected 0.05 < z < 0.20 
galaxies from SDSS DR4 with |Kauffmann et al.||2003] mass esti- 
mates. 

galaxies' stellar mass, although the contribution of the 
asymptotic giant branch remains difficult to model. We 
found, however, that adding 2MASS JHK magnitudes 
to our fits made only small changes to our stellar mass 
estimates, and we do not include them in our stellar mass 
fits. UV observations may be helpful in the cases of 
galaxies with both young and old stellar populations to 
constrain recent star formation history, and we incorpo- 
rate into our an alysis the available stellar masses from 
INeill et al.ll2009l that were fit to GALEX UV as well as 
S DSS u'g'r'i'z' magn itudes. 

Sulli van et al.|2006] carefully studied the effects of more 
complicated star formation histories as well as uncertain 
dust extinction on stellar masses estimates calculated us- 
ing a similar set of PEGASE2 SEDs. They conclude that 
the effect of model uncertainties on stellar masses is small 
althoug h st ar formation rates are more strongly affected. 

Table [12] lists our galaxy measurements. 

5.2.2. Comparison with Kauffmann Masses 

For 10000 randomly selected SDSS DR4 galaxies, we 
compare stellar mass estimates calculated from PE- 
GASE2 fits to u'g'r'i'z' photometry against sp ectro- 
scopic mass estimates from|Kauffmann et al.|20 03 , which 



rely on the |Bruzual fe Charlot||2003 



stellar population 



synthesis model. [Kauffmann et ai. 2003 measure the 
spectral indices D„(4000), the 4000A break, and B5 A , 
stellar Balmer absorption, to constrain star formation 
history. With a reliable star formation history, internal 
dust extinction can be estimated by comparing the model 
prediction to the broadband photometric colors. From 



10 300 300 Myr spheroidal 

2 100 500 Myr spheroidal 

0.66 500 inclination- averaged 

0.4 1000 inclination-averaged 

0.2 1000 inclination-averaged 

0.1 2000 inclination- averaged 

Note. - Summary of PEGASE2 scenarios. 
SFR=^x WLgas where v has units Gyr~ . M gas 
is the gas density. Infall time scales in Myrs. 
Spheroidal dust extinction is fitted with a King's 
profile where the dust density is a power of the 
stellar density. For inclination-averaged extinction, 
dust is placed throughout a plane-parallel slab. 

star formation histories consistent with the spectral in- 
dices and a n estimate of internal extinction, iKauffmann 
|et al.|2003] find each galaxy's mass-to- light ratio and cal- 
culate a mass from the z'-band magnitude. 

Fig ure [T] shows a comparison between Kauffmann et al. 
2003 masses and our mass estimates for 10000 0.05 < 
z < 0.2 galaxies randomly selected from the SDSS DR4 
fit from petroMag catalog magnitudes. Galaxies beyond 
z = 0.05 are sufficiently small and faint enough that the 
SDSS pipeline photometry is reliable, and the SDSS 3" 
fiber samples a significant fraction of the galaxy light. 
The plot shows that the two mass estimates are consis- 
tent, with a root-mean-squared dispersion of 0.15 dex 
and a mean bias of 0.033 dex. 



6. RESULTS 

6.1. Comparison of Low Redshift SN Host Masses to 
SNLS lst-Year Host Masses 

Figure [2] compares host masses of SN in our low red- 
shift sample against those of SN discovered by the Su- 
pernova Legacy Survey (SNL S) in its first year as pub- 
lished by |Howell et at]|2009l The host masses of SNLS 
SN la plott ed in Figure l2| were calculated by |Sullivan 
et al.|[2006| using PEGASE2 stellar population synthe- 
sis models with a somewhat different set of star forma- 
tion histories than we use to analyze our sample of low- 
redshift host galaxies. Our choice of star formation his- 
tories have a 0.15 dex RMS dispersion in residuals from 
Kauffman n et al.|200 3 spectroscopic estimates, while the 
SNLS masses have a less precise RMS dispersion of 0.17 
dex. Nonetheless, SNLS and our mass estimates agree in 
the mean: after fitting stellar masses to u'g'r'i'z' pho- 
tometry of 10000 randomly selected DR4 galaxies with 
our set of templates and the SNLS templates, we find a 
mean offset of 0.009 dex and an RMS of 0.12 dex. Both 
SNLS and this paper model host galaxy ugriz magni- 
tudes, although these correspond to increasingly blue re- 
gions of rest-frame galaxy SEDs for more distant SNLS 
hosts. 

Compared to the SNLS sample, our low redshift host 
galaxy sample has a much higher percentage of hig h mass 
galaxies. This result agrees with |Neill et al.|2009[ who fit 
many of the gala xies in our sample with the PEGASE2 
templates used by |Howell et al.|2009| The high fraction of 
massive galaxies in our sample is likely largely explained 
by the fact that many of our 0.015 < z < 0.08 SN were 
discovered by search programs that regularly target lu- 
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miaous, massive local galaxies. The Lick Observatory 
Supernova Search (LOSS), for instance, cycles through a 
set of host galaxy targets every 3-4 days ( Filippenko et al 



20 



2001 1. Of the 70 galaxies in our sample, 17 were discov- 
ered by the Lick Observatory Supernova Search (LOSS), 
11 were discovered by the Puckett Observatory, and 7 
were discovered by the Lick Observatory and Tenegra 
Observatory Supernova Search (LOTOSS), all of which 
perform targeted SN searches. The SNLS, in contrast, is 
an untargeted survey which observed four fields during 
dark and gray t ime every 3-5 nights for 5-6 lunar phase 



cycles per year (Balland et al. 2009]). Another effect fa- 



voring the discovery of SJN la in massive galaxies at low 
redshift is the fact that most SN la in massive elliptical 
galaxies are fast-declining. These faint SN are less likely 
to be discovered at higher redshifts. 

Given the correlation between the star formation rate 
and the SN la rate, we expect that the SN la rate will 
increase at higher redshift where the average star forma- 
tion rate is higher. Current understanding of star forma- 
tion histories suggests that more massive galaxies form 
their stars earlier and over a shorter period of ti me, an 



effect known as 'downsizing' (Cowie et al. 1996). This 



suggests that the rate of SN la with shorter delay times 
in more massive galaxies may be elevated at higher red- 
shifts. The relatively constant number density of galax- 
ies w ith masses greater th an 10 10 8 M Q since z = 0.75 
(i.e. jPozzetti et al. 2009) suggests that an intermedi- 
ate redshift survey may discover a greater fraction of 
slow-declining SN la in galaxies with masses greater than 
10 108 Mq compared to an untargeted low redshift sur- 
vey. 



6.2. Hubble Residuals and Errors 

In this paper, we define a Hubble residual as 
HR = fj^sN — (J-z so that a positive residual indicates that 
a SN la is fainter after light curve correction than ex- 
pected for the bes t-fit cosmology. As calculated by 
|Hicken et al.|2009b[ the error in HR is the sum in quadra- 
ture of errors in hsn and [i z . Error in /isn includes un- 
certainties from the light curve photometry, extinction, 
the light curve fitter's model coefficients, and an 'intrin- 
sic uncertainty' that yields xt ~ 1 f° r cosmology fits, 
while the uncertainty in [i z comes from peculiar velocity 
(with a 400 km s _1 dispersio n). Here our sampl e is a 
subset of the SN that pass the Hicken et al.|2009b] 'mini- 
m al cuts,' and we remov e the 'intrinsic uncertainty' used 
by |Hicken et al . 2009b, We add uncertainties that give 
xt ~ 1 for the subsamples corresponding to each fitter 
(0.14 mag for MLCS17, 0.11 mag for MLCS31, 0.13 mag 
for SALT, and 0.13 mag for SALT2). 

We use t he 400 km s -1 pecul iar velocity dispersion 
adopted by |Hicken et al.|[2009b| to calculate errors on 
galaxy properties and the measurement error covariances 
between Hubble residuals and galaxy properties. A pecu- 
liar velocity towards us makes a SN appear fainter than 
expected for its redshift and also leads to an underesti- 
mation of the host galaxy's stellar mass and size. Con- 
versely, a peculiar velocity away from us makes a SN 
appear brighter than expected for its redshift and also 
leads to an overestimation of the host galaxy's stellar 
mass and size, producing measurement error covariances 
between the SN Hubble residual and host properties. 
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Fig. 2. — Stellar mass comparison between hosts from this paper's 
sample and from the SNLS lst-year survey. Both sets of masses 
were derived by marginalizing over fits of PEGASE2 stellar popu- 
lation synthesis models to ugriz galaxy photometry, although we 
use a different set of star formation histories that produces better 
agreement with spectro scopically-dctcrmined mass estimates from 
IKauffmamT et al. 2003] These two sets of PEGASE2 templates 
produce mass estimates with a mean residual offset of 0.009 dex 
and a root-mean-square spread of 0.12 dex. 

6.3. Accounting for Calibration Error that Correlates 
with Light Curve Parameters 

We want to identify trends in Hubble residuals with 
host galaxy size or mass that cannot be equally well ac- 
counted for by a light curve-dependent calibration error. 
Although host properties and light curve parameters are 



correlated (i.e. Hamuy et al. 19961, their correlation is 



not perfect. The presence of SN la with a diversity of 
light curves in host galaxies within a given size or mass 
range allows us to disentangle trends in Hubble residuals 
with light curve parameters from trends with host galaxy 
properties. We fit for linear trends with a host galaxy 
property (size or mass) and two light curve parameters 
simultaneously so that, for example, 

HR = a x (host mass) + /? x (stretch) + 7 x (color) +5 (3) 

in the case of the SALT and SALT2 fitters where HR 
is the Hubble residual. If the data cannot differenti- 
ate between a Hubble residuals trend with host mass 
or stretch, for instance, then only the linear combination 
of a and (3 will be constrained. Consequently, a and 
(3 can cancel each other out, and they will individually 
be highly uncertain. If instead a fit finds strong con- 
straints on a and (3, then any trends in host mass and 
stretch are independent systematic effects in the data. 
Because host properties and light curve parameters are 
correlated, we marginalize over uncertainties in the co- 
efficients of trends with light curve parameters when we 
report uncertainties in the coefficients trends with host 
size or mass. This marginalization occurs during cal- 
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Fig. 3.— Trends in Hubble residuals with MLCS17 A and SALT2 xi. Of our MCMC realizations, 0.7% (2.4<r; 60) have slopes greater 
than zero for MLCS17 and 2.3% (2. 0<r; 58) have slopes less than zero for SALT2. In parentheses is the significance of a non-zero slope 
followed by the number of SN included in the fit. In our analysis of possible trends in Hubble residuals with host properties, we calculate 
statistics where we marginalize over linear trends with A and stretch. Trends in Hubble residuals with Ay or c trend (7) for MLCS17 and 
SALT2 have slopes is consistent with zero. Of our MCMC realizations, 36.4% (0.4cr; 60) have a slope with Ay less than zero for MLCS17 
and 13.7% (1.1c; 58) have a slope with c less than zero for SALT2. 



culation of posterior probabilities in our Markov Chain 
Monte Carlo (MC MC) analys is implemented with the 
LINMIX package (|Kelly||2067|). Our fits take into ac- 



count measurement errors as well as measurement error 
covariance that arises from peculiar velocity. 

Here we plot trends in Hubble residuals with light 
curve parameters to show the reader their size and di- 
rection. For brevity, we plot only the trends in MLCS17 
Hubble residuals with A and SALT2 residuals with x\ 
in Figure [3] Plots of Hubble residuals versus Ay for 
MLCS17 and c for SALT2 are not shown because best- 
fit trends have slopes consistent with zero. The t rend in 
Hubble resi duals with MLCS17 A, as discussed by |Hicken| 
|et al.|[2009a| may reflect the negative quadratic term in 
A needed to accurately calibrate the quickly decaying, 
1991bg-like SN la. Removing the 1991bg-like SN from 
the training set may result in a positive quadratic term 
which provides a bette r fit as a function of A for 'normal' 
SN la [see Figure 26 of |Hicken et al |2009a| . Table [3] lists 
statistics for all four light curve fitters. 

6.4. Hubble Residuals and Hosts' Physical Sizes 

We find a trend in Hubble residuals with the ra- 
dius of the circular aperture enclosing 90% of the 
galaxy i-band light, calculated from SExtractor parame- 
ter FLUX_RADIUS with PHOT_FLUXFRAC=0.9. Fig- 
ure |] plots the 90% host radii against MLCS17 and 
SALT2 Hubble residuals. The outlying data point with 
large host radius in the panels is SN 2003ic whose host 
is the cD galaxy at the center of the galaxy cluster Abell 
85. 

First, we fit only for a trend in Hubble residuals with 
host galaxy size (a) while holding ^=7=0. We fit a Gaus- 
sian to the posterior probability distribution and use the 
standard deviation to calculate the a significance of a 
non-zero slope. The significance of a non-zero slope (a) 



ranges from 2.6-2.9(7 depending on the light curve fitter. 
Full statistics are summarized in Tables |4j [5j [7j and [8j 
Separating the SN by host galaxy 90% radius into two 
equals-sized bins, the larger-host bin has more negative 
Hubble residuals according to the bins' weighted aver- 
ages. The significance of a difference between the bins' 
weighted averages ranges from 2.3-3.2(7 among the light 
curve fitters. 

We then simultaneously fit for trends in Hubble resid- 
uals with host size (a) and light curve parameters (/3,7) 
A and Ay for MLCS2k2 fitters and stretch and color 
for SALT and SALT2. Marginalizing over /3 and 7, the 
significance of a non-zero slope (a) ranges from 1.7-2.4cr 
depending on the light curve fitter (2.0-2.3<7 excluding 
MLCS31). The joint posterior probability distributions 
do not reveal strong degeneracies between the fit coeffi- 
cients. Removing the best-fitting trends with light curve 
parameters (fi-Ti) from the Hubble residuals after the 
simultaneous fit, a 1.9-2.4cr difference between the bin 
weighted averages persists. 

The cD host galaxy of SN 2003ic is significantly larger 
than other hosts, and, to test the effect of this data point 
on our statistics, we repeat the analysis after removing 
it from the sample. The bin weighted average differences 
have significances of 1.9-3.1(7, although a linear fit to the 
data with a negative slope has only 1.4-2.2er significance 
(2.0-2.2cr without SALT). 

Tables [6] and [9] show that improvement in the standard 
deviation of Hubble residuals after removing the host- 
dependent trend is approximately 5-10%. 

6.5. Hubble Residuals and Galaxy Masses 

Mass is a fundamental property of galaxies which helps 
determine their star formation histories and chemical en- 
richment. Spectroscopic metallicity measurements have 
shown that stellar mass correlates with galaxy gas-phase 
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Fig. 4. — Trends in MLCS17 and SALT2 Hubble residual with host radius enclosing 90% of i-band light. SN host radius correlates with 
Hubble residual for both (a) MLCS17 and (b) SALT2, with the SN that occur in larger galaxies having significantly more negative Hubble 
residuals. The weighted averages of two bins in host size have a magnitude difference of 0.11 (2.3<r; 60) for MLCS17 and 0.11 (2.6<j; 62) for 
SALT2. When we fit for the trend with host radius (a) while holding /3=7=0, only 0.5% (2.6cr; 60) of slopes drawn from an MCMC 
analysis have greater than zero slope for MLCS17 and 0.4% (2.7<r; 58) for SALT2. In parentheses is the significance of a non-zero slope 
followed by the number of SN included in the fit. Posterior slope distributions are plotted in the upper panels. The outlying data point is 
SN 2003ic, hosted by the cD galaxy at the center of the Abell 85 galaxy cluster. With SN 2003ic removed, the bin weighted averages differ 
in magnitude by 0.12 (2.7a) for MLCS17 and 0.11 (2.6<r) for SALT2 while 1.6% (2.2ct; 59) of slopes drawn from an MCMC analysis are 
greater than zero for MLCS17 and 2.2% (2.0cr; 61) for SALT2. 
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Fig. 5. — Trends in MLCS17 and SALT2 Hubble residual with host stellar mass (log M> 9.5). SN with more massive hosts have Hubble 
residuals that are more negative in the cases of (a) MLCS17 and (b) SALT2, consistent with the trend evident in Figure H] Host stellar 
masses were measured using fits to ugriz photometry with the PEGASE2 stellar population synthesis templates. The weighted averages 
of two bins separated by their masses yield a magnitude difference of 0.11 (2.3<r; 58) for MLCS17 and 0.11 (2.5ct; 60) for SALT2. When we 
fit for the trend with host stellar mass (a) while holding /3=7=0, only 0.2% (2.8<r; 58) of slopes drawn from a MCMC analysis are greater 
than zero for MLCS17 and 0.4% (3. Oct; 60) for SALT2. In parentheses is the significance of a non-zero slope followed by the number of SN 
included in the fit. The upper panels plot the posterior slope distributions. 
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( |Tremonti et al.||2004[ ) and stellar ( |Gallazzi et al.||2005 l 
metal content, reflecting the fact that it is easier tor met- 
als to escape from the shallower gravitational potentials 
of less massive galaxies. As seen in Figure [2j only 2 of 
our SN host galaxie s have masses less than 10 M . 
|Timmes et al.|[2003| argue that the effect of metallicity 
on SN la luminosity depends non-linearly on the pro- 
genitor's metallicity, which might mean that any effect 
on Hubble residuals would only become strong for SN 
la in massive galaxies. We restrict our sample for our 
primary analysis to SN with galaxy hosts more massive 
than 10 9 - 5 M where we have a sufficient number of SN 
with massive hosts to constrain any systematic trend. 

We plot galaxy stellar mass against Hubble residuals in 
Figure [5] for the MLCS17 and SALT2 light curve fitters, 
excludingthe 2 SN with the least massive hosts. Tables|4j 
[5j[7j andJS] summarize all of our statistical tests. First, we 
nt only for a trend in Hubble residuals with host galaxy 
mass (a) while holding /3=7=0. The significance of a 
non-zero slope (a) ranges from 2.3-2.8cr depending on 
the light curve fitter. Separating the SN by host galaxy 
mass into two equal-sized bins, the more massive host 
bin has more negative Hubble residuals according to the 
bins' weighted averages. The significance of a difference 
between the bins' weighted averages ranges from 2.2-2.8cr 
among the light curve fitters. 

We then simultaneously fit for trends in Hubble resid- 
uals with host mass (a) and light curve parameters (/3,7) 
A and Ay for MLCS2k2 fitters and stretch and color for 
SALT and SALT2. Marginalizing over (3 and 7, the sig- 
nificance of a non-zero slope (a) ranges from 2.0-2.6cr 
depending on the light curve fitter. The joint poste- 
rior probability distributions do not reveal strong degen- 
eracies between the fit coefficients. Removing the best- 
fitting trends with light curve parameters (/3,7) from the 
Hubble residuals after the simultaneous fit, a 1.6-2.6cr 
difference between the bin weighted averages persists. 

Tables H] and [9] show that the reduction in the standard 
deviation of Hubble residuals after removing the host- 
dependent trend is approximately 5-10%. 

6.5.1. Weill et al.\[20M Masses 

We now repeat the anal yses using stellar mass esti- 
mates from Neill et al.|2009| which were fit using GALEX 
UV measurements in a ddition to SPSS u 'g'r'i'z' mag- 
nitudes. Because each Neill et al. 2009 host galaxy is 



required to have UV as well as optical measurements, 
|Neill et al.|[2009| measured masses for only 49 of the 70 
SN in our sample with MLCS17 A v < 0.5. Fitting only 
for a trend in Hubble residuals with host galaxy mass 
(a) while holding /3=7=0, the significance of a non-zero 
slope (a) ranges from 2.3-2.7er depending on the light 
curve fitter. The significance of a difference between the 
bins' weighted averages ranges from 1.8-3.2cr among the 
light curve fitters. 

Simultaneously fitting for and marginalizing over linear 
trends (/3,7) with light curve parameters, the significance 
of a non-zero slope (a) is 1.9<r for all light curve fitters. 
The joint posterior probability distributions do not reveal 
strong degeneracies between the fit coefficients. Remov- 
ing the best-fitting trends with light curve parameters 
(f3,j) from the Hubble residuals after the simultaneous 
fit, a 1.6-2. 5er (2. 1-2. 5a without SALT2) difference be- 
tween the bin weighted averages remains. 



Tables [6] and [9] show the improvement in the standard 
deviation of Hubble residuals after removing the host- 
dependent trend. This results in improvements of up to 
15% in standard deviation. 

6.6. Information Criteria 

The Bayesian information criterion (BIC) and the 
Akaikc information criterion (AIC) help to evaluate 
whether an additional model parameter can be justified 
by the improvement it produces in the % 2 goodness-of- 
fit statistic. BIC derives from an approximation to the 
Bayes factor, the ratio of the probabilities of two models 
given equal prior probability. For any two prospective 
models, the model with the lower value of BI C is the 
prefe rred one. With Gaussian errors, the BIC (Schwarz 
1978 1 is given as, 



BIC = x 2 + fclnJV 



(4) 



where k is the number of parameters and N is 
the number of data points in the fit. A compari- 
son can be made between two mo dels by comp uting 
ABIC = Ax 2 + AklnN. The AIC ( | Akaike| [l974| ) uses 
Kullback-Leibler information entropy as a basis to com- 
pare different models of a given data set. The AIC e is a 



versio n of the AIC corrected for small data sets ( Sugiura 
T978|, 



AIC C = x +2k 



2k(k + l) 
N-k-1 



(5) 



where k is the number of parameters and N is the number 
of data points used in fit. A decrease of 2 in the infor- 
mation criterion is positive evidence against the model 
with higher information criterion, while a decrease of 6 
is strong p ositive evidence against the model with large r 
value (e.g. |Kass & Raftery|1995j |Mukherjee et al.|1998[ ). 
To calculate these statistics, we first add an 'intrinsic un- 
certainty' in quadrature to the Hubble residual error on 
each SN so that xt — 1 f° r a fit with no trends. We 
calculate the BIC and AIC C for three separate models: 
no trend (S), a trend with light curve shape parameter 
(f3,S), and trends with light curve shape parameter and 
host galaxy property (a,/3,<5). Only differe nces in in- 
formation criteria are significant, and Tables |l0| and 11 
list the change in the information criteria relative to the 
model with no trend (5). An Ay/c trend (7) is not in- 
cluded because its slope is consistent with zero and its 
inclusion results in no significant reduction in the fit \ 2 - 
We are interested in the changes in the BIC and AIC C 
after adding a host dependent trend (a) to a model that 
already includes a trend with light curve shape (/3,<5). 
For all light curve fitters, an added trend with host size 
produces ABIC < —2, providing evidence in favor of 
the additional trend. For trends with host stellar mass, 
ABIC < -2 exc ept for SALT. ML CS17 and SALT2 have 
ABIC < -2 for (Neill et aLl[2009l stellar masses. Only if 



ABIC > 2 would there be evidence against adding a host 
dependent term, and in no case is ABIC > 2. AAIC C is 
more negative in all cases than ABIC because BIC has 
a stronger penalty for each additional parameter. In the 
cases of all trends and all light curve fitters, AAIC C < —2 
and largely is less than -4, providing moderate to strong 
positive evidence for a host dependent term. 
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Fig. 6. — Hubble residual trend persists after i'-band surface brightness cut. (a) plots the asinh magnitude %' surface brightness in a 0.2 
kpc aperture at the SN locations. We make a cut at a surface brightness of 21 asinh mag arcsecond" 1 indicated by the dashed horizontal 
line. The surface brightness distribution of SN locations fainter than this threshold is similar across host masses, (b) shows that the 
correlation between Hubble residual and host mass exists for these SN that occur in lower surface brightness locations, indicating that a 
Malmquist-like effect is not likely to be responsible for the observed trends in Hubble residuals with host galaxy mass. 



6.7. Possible Systematic Biases 
6.7.1. Malmquist Bias 

Photographic searches had a well-known bias against 
SN detection near galaxy centers where images were of- 
ten saturated ( Shaw|1979 1 . The much improved dynamic 
range of CCD cameras used by modern SN search pro- 
grams and digital image subtraction enables detection 
in even high surface-brightness regions. More massive 
galaxies on average have regions of higher surface bright- 
ness, and we investigate whether a detection bias could 
explain the ~0.10 mag effect associated with the galaxy 
mass. To account for a ~0.10 mag trend, SN searches 
would have to miss a high percentage of only slightly 
fainter explosions in brighter, more massive hosts. In 
particular, the detection efficiency in massive hosts would 
have to decrease to near zero over a ~10% range in SN 
brightness. However, if there were such a steep drop 
in efficiency for brighter SN with intermediate decline 
rates, SN searches would detect few if any fast-declining 
SN la in massive galaxies because their peak intrinsic 
brightness can be at least 50% fainter than the SN with 
intermediate decline rates detected in these galaxies. Be- 
cause we do see SN la with both faster and slower decline 
rates in massive galaxies, it seems unlikely that a detec- 
tion bias explains the effect with host stellar mass we 
observe. 

To test for the possible effect of Malmquist bias, we ap- 
ply a cut to remove SN la detected in high-surface bright- 
ness regions so that the SN in low and high mass hosts 
have similar distributions of galaxy surface brightness at 
the locations of their explosions. Figure |6][a) plots the 
j'-band surface brightness in a 0.2 kpc aperture at the 
SN la location. Applying a cut at a surface brightness of 
21 asinh mag arcsec^ 1 , shown by the dashed horizontal 
line, leaves surface brightness distributions that are sim- 
ilar for low and high mass hosts. Figure [6jb) shows that 
a trend in Hubble residuals with host mass persists for 



this subset, suggesting that Malmquist bias is likely not 
responsible for the trend we observe. 

6.7.2. Regional Velocity Flows 

Peculiar galaxy velocities arising from departures from 
the Hubble flow could have a systematic effect on Hubble 
residuals. We perform the simple test of splitting our SN 
sample into lower and higher redshift groups. Both of 
these samples are consistent with size and mass depen- 
dent trends. 

6.7.3. Potential Correlations Analyzed 

The initial objective of our analysis was to constrain 
any correlation between the host galaxy stellar masses 
and Hubble residuals o f our nearby SN sa mple to pro- 
vide comparison to the |Howell et al.|[2009| SNLS study. 
We subsequently expanded our search to include poten- 
tial correlations of Hubble residuals with additional host 
properties including galaxy color as well as the location 
of the SN in the host galaxy. Because galaxy properties 
including stellar mass, size, and color correlate with each 
other, the number of independent parameter searches is 
difficult to estimate precisely but likely is ~2-3. The 
probability of finding at least one correlation with sig- 
nificance p in an analysis of n potential correlations is 
p' = 1 — (1 — p) n . Accordingly, a search for correlations 
with 2 independent variables reduces the significance of 
a 2.5a result to ^2.25cr, while a search for correlations 
with 3 independent variables reduces the significance of 
a 2.5<7 result to ~2.1<7. 

7. DISCUSSION 

Given the correlation between light curve parameters 
and host galaxy properties, we have been careful to in- 
vestigate whether the trends we detect with host size and 
mass could instead be accounted for by trends with SN 
la light curve shape or color (or dust). SN with a range 
of light curve properties are present in the high and low 
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mass ends of the host galaxy distribution, and this allows 
us to disentangle potential effects of size, mass, and light 
curve parameters. By simultaneously fitting for and then 
marginalizing over trends with light curve parameters in 
our fits for host-dependent trends, we have shown that 
calibration errors that correlate with light curve prop- 
erties do not account for the trends we find in Hubble 
residuals with host galaxy size and mass. 

Host galaxy metallicity, stellar age, and extinction by 
dust correlate with galaxy mass and may be factors re- 
sponsible for the variation we observe in the SN la width- 
color-luminosity relation. According to the |Timmes et ah] 
|2003| model, the effect of metallicity on SN la lumi- 
nosities depends non-linear ly on metallicity and would 
be strongest at the hi gh metallicities of the massive 
galaxies in our sam ple. Kasen et al.||2009| included the 
ITimmes et al. 120031 effect in SN la simulations and found 
that increased metallicity would alter the SN la width- 
luminosity relation so that SN la with intermediate to 
slow decline rates would be intrinsically brighter after 
light curve correction. Such a correction would be con- 
sistent with the direction of the trend we see in Hubble 
residuals with increasing host masses, at l east for SN 
with slower decay rates. |Kasen et al.||2009| did not ac- 
count for the possible effects of metals on the dynamics 
of the explosion, or on the white dwarf's structure. We 
note that an interesting feature of Figure [4] is the appar- 
ent lack of a strong trend in the Hubble residuals with 
radius for SN la in host galaxies with sizes greater than 
15 kpc. 

That the hosts of our low redshift SN la are substan- 
tially more massive than hosts of 0.2 < z < 0.75 SNLS 
SN helps to interpret existing studies of Hubble resid- 
uals and host gala xy properties. Figure [2] shows that 
|Howell et aT] [20091 analyzing SNLS SN, fit for a linear 
trend in Hubble residuals over a much larger range of host 
masses than we do and have comparatively few galaxies 
with masses greater than 10 11 M©, perhaps helping to 
explain the fact that they found no strong trend in Hub- 
ble residuals with host metallicity. Gallagher et al.|2005 
studied SN la in star-forming and therefore likely less 
massive low redshift galaxies, finding no significant trend 
in Hu bble residuals with host metallicity. Gallagher et al. 



2008 however, studied low redshift E/S0 galaxies, which 
are likely similar to the more massive hosts in our sample 
(four hosts overlap with our sample), and found a 98% 

confide nce in a trend in host me tallicity. 

Both [Gallagher et al |2008|and|Howeil et al.|2009| favor 



trends (although the Howell et al.||2009| trend only has 
l-o" significance) where galaxies with higher metallicities 
host SN that are brighter than expected. If the more 
massive galaxies in our sample have higher metallicities, 
the direction of the trends we find agrees with the direc- 
tion of the trend in Hubble residuals with host metallic- 
ity found by |Gallagher"et~aLl|2008l The |Neill et aTl2009] 
trend in Hubble residuals with host age exists only in a 
subset of SN whose hosts have low extinction but whose 
light curve colors are not necessarily consistent with low 
extinction, making comparison with our trend difficult. 

While we could in theory use an empirical mass- 
metallicity relation to test di rectly for the effect of metal- 
licity in our sample following |Howell et al.|2009| and [Neill| 
|et al.|[2009[ it is not clear which mass-metallicity rela- 
tion is appropriate. The fact that SN la progenitors are 



likely to be younger in later-type hosts may mean that 
their metallicities are similar to the gas-phase metallic- 
ity while SN la progenitors in older, early-type galaxies 
may track the stellar metallicity. The slope of mass- 
metallicity relation also becomes flatter for galaxies with 
approximately 10 10 5 M , making mass a less sensitive 
proxy for metallicity. 

7.1. Strong Cosmology Implications 

Low redshift SN are essential to accurate constraints 
on the dark energy equation of state parameter, w, be- 
cause they anchor the redshift-distance relation and set 
the expectation for the brightnesses of high redshift SN 
according to different cosmologies. The discovery that 
the universe is accelerating in its expansion originated 
from the observation that high redshift SN la were fainter 
than expected for a matter- only universe ( Riess et al 



1998 Pcrlmutter et al. 



1999). Low redshift SJN are also 



used to train the MLCS2k2 light fitters and consequently 
influence distance measurements to high redshift SN. As 
Figure [2] demonstrates, only a small percentage of SN dis- 
covereaby the SNLS flux-limited survey have host galaxy 
mass greater than 10 10 " 8 M . However, half of the low 
redshift sample has hosts more massive than 10 10 8 M , 
and this paper has shown that these SN have Hubble 
residuals that are on average more negative by ~0.10 
mag. 

To examine the sensitivity of the best-fitting value 
of w to the host-dependent calibration error, we divide 
our sample into two sets according to whether the SN 
host has a mass greater or less than 10 108 M Q . Us- 
ing the wfit c osmology fitter from t he Supernova Anal- 
ysis package (| Kessler et al . 2009b) and incorpor ating 
baryon acoustic oscillation (Eiscnst ein et al. 2005) and 
cosmic microwave background ( |Komatsu et al. 2009 ) 
measurements, the best-fitting values of w differ. When 
fit in combination with 180 higher-redshift ESSENCE, 
SNLS, and HigherZ SN, the 30 SN with lower mass 
hosts yield 1 + w = 0.22;°;^g while the 30 nearby SN 
with higher mass hosts yield 1 + w = — 0. 03^^43 from 
MLCS17 measurements. The difference between the two 
estimates of w is likely more significant than indicated by 
the overlap between their constraints, because they in- 
clude 180 higher-redshift SN in common. The 30 nearby 
hosts with masses less than 10 108 M Q in our sample are 
more massive than a majority of SNLS host galaxies, so 
we note that we have not constructed a sample of SN 
la with hosts of similar masses across all redshifts. The 
host-dependent shift in w is greater than the statistical 
error on w and implies that host galaxy or progenitor- 
dependent effects have an unappreciated and strong in- 
fluence on even current SN cosmological estimates. 

8. CONCLUSIONS AND SUMMARY 

The hosts of our sample of 70 low redshift (0.015 < 
z < 0.08) SN are substantially more massive than hosts 
of the intermediate redshift (0.2 < z < 0.75) SNLS SN. 
In our sample, SN la that occur in larger, more massive 
galaxies are ~10% brighter than SN la found in smaller, 
less massive galaxies after light curve correction when we 
calculate distance moduli using the MLCS17, MLCS31, 
SALT and SALT2 fitters. When we fit simultaneously 
for trends with light curve parameters in addition to 
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host galaxy properties, we find that these host-dependent 
trends in Hubble residuals cannot be accounted for by 
light curve-dependent calibration error. 

The variation in the SN la width-color-luminosity rela- 
tion with host size and mass implied by these trends may 
reflect the correlation between galaxy mass and metallic- 
ity, increased average ages of progenitors in more massive 
hosts, or another environmental difference such as dust 
extinction. Sensitivity of SN la luminosities to progenitor 
metallicity may become greatest at the high me tallicities 
found in th e most massive hosts in our sample (Timmes 
etaL |200^ . 



IN ear by SN la detected by targeted searches should be 
treated carefully when included in cosmology fits because 
their host galaxies are generally substantially more mas- 
sive than those of SN detected by flux-limited surveys at 
higher redshifts. The 30 nearby SN with host masses less 



than 10 M in our sample yield 1 



w 



0.22 



+0.152 
^0.108 



while the set of 30 nearby SN with more massive hosts 



yield 1 + w = -0.03 



+0.217 
-0.143 



with MLCS17 distances when 



fit in combination with 180 higher-redshift ESSENCE, 
SNLS, and HigherZ SN. By combining information from 
host galaxy measurements with light curve fits, distances 
to SN la may be more accurately estimated and used to 
measure cosmological parameters with better control of 
systematics. 
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TABLE 3 

Analysis of Potential Light Curve Dependence of Hubble Residuals 



MLCS17 


MLCS31 


SALT 


SALT2 


A,s,xi 0.7% (2.4cr; 60) 
A v ,c 36.4% (0.4o-; 60) 


3.6% (1.9a; 58) 
0.4% (2.4cr; 58) 


13.8% (l.ler; 26) 
20.3% (0.8cr; 48) 


2.3% (2.0a; 58) 
13.7% (l.lcr; 58) 



Note. — Each statistic is the percentage of MCMC realizations that have slopes with 
opposite sign from the slope with greatest posterior probability. In parentheses is the 
significance of a non-zero slope followed by the number of SN included in the fit. 



TABLE 4 

Magnitude Difference Between Bin Weighted Averages of Hubble Residuals 



MLCS17 MLCS31 SALT SALT2 



90% Radius 0.11 (2.3er; 60) 0.17 (3.2cr; 58) 0.087 (1.9er; 53) 0.11 (2.6a; 62) 
90% Radius (<30 kpc) 0.12 (2.7a) 0.17 (3.1a) 0.087 (1.9er) 0.11 (2.6a) 

log M (>9.5) 0.11 (2.3a; 58) 0.12 (2.2<r; 56) 0.13 (2.8cr; 51) 0.11 (2.5cr; 60) 
log M 0.10 (2.2a) 0.13 (2Aa) 0.11 (2.5a) 0.051 (I.Oct) 

log M (>9.5; GALEX, Neill 09) 0.17 (3.2cr; 45) 0.19 (3.0a; 45) 0.14 (2.7cr; 40) 0.094 (1.8<r; 47) 
log M (GALEX, Neill 09) ' 0.19 (4.0ct) 0.21 (3Aa) 0.14 (2.7a) 0.11 (2.2a) 

Note. — We separate the SN into two equal-sized bins according to host property and calculate the magnitude 
difference between the weighted averages of Hubble residuals of SN in each bin. In parentheses is the significance 
followed by the number of SN included in the fit. 



TABLE 5 

Significance of Non-Zero Slope in Hubble Residuals 



MLCS17 



MLCS31 



SALT 



SALT2 



90% Radius 
90% Radius (<30 kpc) 
log M (>9.5) 
log M 

log M (>9.5; GALEX, Neill 09) 
log M (GALEX, Neill 09) 



0.5% (2.6a; 60) 
1.6% (2.2a; 59) 
0.2% (2.8ct; 58) 
0.4% (2.4cr; 60) 
0.4% (2.8cr; 45) 
0.01% (3.9<r; 46) 



0.4% (2.7a; 58) 
1.4% (2.2o-; 57) 
0.2% (2.7a; 56) 
0.7% (2.4o-; 58) 
0.7% (2Aa; 45) 
0.04% (3.0<r; 46) 



0.2% (2.7a; 53) 
8.0% (1.4a; 52) 
1.1% (2.3ct; 51) 
29.9% (0.5a; 53) 
0.3% (2.7o-; 40) 
0.7% (2.5a; 40) 



0.2% (2.8cr; 62) 
2.2% (2.0ct; 61) 
0.4% (3.0<t; 60) 
2.3% (2.0a; 62) 
0.7% (2.6cr; 47) 
0.03% (3.5<x; 48) 



Note. — Significance of a non-zero slope. Here we fit for the trend with host radius (a) while holding /3=7=0. 
Each statistic is the percentage of MCMC realizations that have slopes with opposite sign from the slope with greatest 
posterior probability. In parentheses is the significance of a non-zero slope followed by the number of SN included in 
the fit. 
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TABLE 6 

Improvement in the Standard Deviation of Hubble Residuals after Subtracting 

Host-Dependent Trend 



MLCS17 



MLCS31 



SALT 



SALT2 



90% Radius 0.193 \ 0.184 0.222 \ 0.210 0.168 \ 0.156 0.185 \ 0.175 

90% Radius (<30 kpc) 0.191 \ 0.185 0.219 \ 0.212 0.161 \ 0.157 0.182 \ 0.177 

log M 0.193 \ 0.182 0.222 \ 0.208 0.168 \ 0.167 0.185 \ 0.178 

log M (>9.5; GALEX, Neill 09) 0.159 \ 0.143 0.192 \ 0.180 0.174 \ 0.159 0.163 \ 0.150 

log M (GALEX, Neill 09) 0.181 \ 0.151 0.213 \ 0.187 0.174 \ 0.158 0.184 \ 0.158 

Note. — First value is the standard deviation in magnitudes before and the second value is the 
standard deviation in magnitudes after subtracting the trend with host property, a X (host property), 
determined from fits where we hold /9=7=0. 



TABLE 7 

Magnitude Difference Between Bin Weighted Averages of Hubble Residuals After Removing Trends 

with Light Curve Parameters 



MLCS17 



MLCS31 



SALT 



SALT2 



90% Radius 
90% Radius (<30 kpc) 
log M (>9.5) 
log M 

log M (>9.5; GALEX, Neill 09) 
log M (GALEX, Neill 09) 



0.089 (1.90-; 60) 
0.10 (2.2a) 

0.093 (1.9a; 58) 
0.083 (1.8er) 

0.13 (2.5cr; 45) 
0.17 (3.2<t) 



0.11 (2.1a; 58) 
0.11 (2.0a) 

0.089 (1.6(j; 56) 
0.085 (1.5a) 

0.13 (2.1a; 45) 
0.17 (2.8a) 



0.082 (1.8<t; 53) 
0.083 (I.80-) 

0.12 (2.6a; 51) 
0.090 (2. Oct) 

0.11 (2.1<t; 40) 
0.11 (2.1a) 



0.10 (2Aa; 62) 
0.10 (2.3a) 

0.10 (2.2cr; 60) 
0.094 (2.1a) 

0.083 (1.6c-; 47) 
0.10 (2.0ct) 



Note. — Here we fit for trends with host property (a) and light curve parameters (/3,7) simultaneously. We 
remove trends with light curve parameters (/3,7) from the Hubble residuals and separate the SN into two equal- 
sized bins according to host property. We then calculate the magnitude difference between the weighted averages 
of Hubble residuals of SN in each bin. In parentheses is the significance followed by the number of SN included in 
the fit. 



TABLE 8 

Significance of Non-Zero Slope in Hubble Residuals After Marginalizing over Light Curve 

Parameters 



MLCS17 



MLCS31 



SALT 



SALT2 



90% Radius 
90% Radius (<30 kpc) 
log M (>9.5) 
log M 

log M (>9.5; GALEX, Neill 09) 
log M (GALEX, Neill 09) 



1.9% (2. 0<t; 60) 

7.0% (1.5cr; 59) 

1.6% (2.2a; 58) 

5.8% (1.6er; 60) 

3.0% (1.9a; 45) 

0.2% (3.0cr; 46) 



2.3% (1.9a; 58) 
7.6% (1.5a; 57) 
3.2% (1.8a; 56) 
11.4% (1.2a; 58) 
3.9% (1.7a; 45) 
0.5% (2.7<x; 46) 



0.6% (2Aa; 53) 
9.0% (1.4ct; 52) 
2.3% (2.0ct; 51) 
38.4% (0.3<r; 53) 
2.0% (2.0a; 40) 
2.1% (2.0ct; 40) 



0.9% (2.3a; 62) 
5.7% (1.6a; 61) 
0.8% (2Aa; 60) 
9.8% (1.2a; 62) 
1.4% (2.1a; 47) 
0.2% (2.8ct; 48) 



Note. — Significance of a non-zero slope. Here we have fit for trends with host property (a) and light curve 
parameters (/3,7) simultaneously. Each statistic is the percentage of MCMC realizations that have slopes with 
opposite sign from the slope with greatest posterior probability. In parentheses is the significance of a non-zero 
slope followed by the number of SN included in the fit. 



TABLE 9 

Improvement in the Standard Deviation of Hubble Residuals after Subtracting 
Host-Dependent Trend from Simultaneous Fit 



MLCS17 



MLCS31 



SALT 



SALT2 



90% Radius 0.193 \ 0.184 0.222 \ 0.210 0.168 \ 0.156 0.185 \ 0.175 

90% Radius (<30 kpc) 0.191 \ 0.185 0.219 \ 0.212 0.161 \ 0.157 0.182 \ 0.177 

log M 0.193 \ 0.183 0.222 \ 0.211 0.168 \ 0.168 0.185 \ 0.180 

log M (>9.5; GALEX, Neill 09) 0.159 \ 0.142 0.192 \ 0.180 0.174 \ 0.158 0.163 \ 0.150 

log M (GALEX, Neill 09) 0.181 \ 0.151 0.213 \ 0.188 0.174 \ 0.158 0.184 \ 0.158 

Note. — First value is the standard deviation in magnitudes before and the second value is the 
standard deviation in magnitudes after subtracting the trend with host property, a X (host property), 
but not the trends with light curve properties (/3,7). Here the trend with host property (a) is determined 
in fits where we simultaneously fit for trends with light curve parameters (/3,7). 
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TABLE 10 
Bayesian Information Criteria 



Parameters 


Host Property 


Light Curve 


MLCS17 


MLCS31 


SALT 


SALT2 


5 






0.0 


0.0 


0.0 


0.0 


0,6 




A/s/xi 


-2.1 


0.6 


-3.2 


-0.3 


a, 13,5 


90% Radius 


A/s/xi 


-7.2 


-3.2 


-7.0 


-4.2 


a, /3,5 


log M (>9.5) 


A/s/xi 


-6.8 


-2.7 


-3.0 


-4.5 


a, /3,5 


log M (>9.5; GALEX, Neill 09) 


A/s/xi 


-5.7 


-0.7 


-3.4 


-2.4 



Note. — BIC increments after adding linear terms (f3,a) to a constant model (<5) of Hubble residuals. Only 
differences between BIC values are meaningful, and models with lower values are favored. A difference between 
and 2 is weak evidence, 2 and 6 evidence, and greater than 6 strong evidence against the model with higher 
BIC. The improvements in \ 2 provide evidence for the addition of trends with host size for all light curve fitter 
and with host mass for MLCS17, MLCS31, and SALT2. An Ay/c trend (7) is not included in this analysis 
because it has a slope consistent with zero and produces no significant reduction in the fit \ 2 ■ 



TABLE 11 

Corrected Akaike Information Criteria 



Parameters 


Host Property 


Light Curve 


MLCS17 


MLCS31 


SALT 


SALT2 


5 


















0,6 




A/s/x-i 


-3.7 


-1.0 


-4.6 


-2.0 


a, 0,5 


90% Radius 


A/s/xi 


-10.3 


-6.3 


-9.8 


-7.4 


a, 0,6 


log M (>9.5) 


A/s/xi 


-9.9 


-5.8 


-5.8 


-7.7 


a, 0,6 1 


og M (>9.5; GALEX, Neill 09) 


A/S/X! 


-8.8 


-3.8 


-6.3 


-5.6 



Note. — AIC C increments after adding linear terms (0,a) to a constant model (5) of Hubble residuals. 
The improvements in \ 2 provide evidence for the addition of trends with host properties (a). AAIC C values 
are generally more negative than ABIC values because BIC includes a stronger penalty for additional fit 
parameters. An Ay/c trend (7) is not included in this analysis because it has a slope consistent with zero and 
produces no significant reduction in the fit x 2 - 



Kelly et al. 



TABLE 12 
Nearby Host Measurements 



SN 


Z CMB 


1)11 '1 r^''lll/l T""l 1 1 1 T 1 c 

V K P L J 


M* — 


M* 


M* + 


1992P 


0.0265 


9.63 


10.31 


10.55 


10.72 


1993ac 


0.0489 


22.32 


11.55 


11.58 


11.62 


1994M 


0.0243 


27.10 


10.99 


11.21 


11.37 


1994S 


0.0152 


10.51 


10.17 


10.68 


10.83 


1994T 


0.0357 


8.15 


10.11 


9.74 


10.74 


1995ac 


0.0488 


17.43 


10.76 


10.99 


11.08 


1996C 


0.0275 


10.02 


9.81 


10.10 


10.30 


1997Y 


0.0166 


8.74 


10.06 


10.59 


10.73 


1997cn 


0.0170 


16.73 


11.52 


11.56 


11.60 


1998ab 


0.0279 


14.58 


10.49 


10.75 


10.93 


1999cc 


0.0315 


21.47 


11.05 


11.23 


11.33 


1999gd 


0.0193 


13.60 


9.89 


10.02 


10.54 


2000cc 


0.0165 


11.79 


10.19 


10.72 


10.85 


2001N 


0.0221 


15.01 


10.81 


10.88 


11.17 


2001V 


0.0162 


16.65 


11.27 


11.30 


11.34 


2001ah 


0.0583 


29.36 


10.86 


11.01 


11.08 


200 Ida 


0.0160 


8.48 


10.10 


10.21 


10.73 


200 lcn 


0.0155 


17.27 


10.45 


11.01 


11.07 


200 lgb 


0.0268 


12.61 


10.42 


10.53 


10.81 


2001ic 


0.0429 


29.91 


11.78 


11.84 


11.88 


2001ic 


0.0312 


15.00 


11.17 


11.26 


11.44 


2002G 


0.0346 


9.76 


10.76 


10.91 


11.03 


2002bf 


0.0249 


10.51 


10.64 


10.84 


11.03 


2002bz 


0.0377 


25.09 


10.98 


11.17 


11.26 


2002ck 


0.0302 


22.08 


10.94 


11.10 


11.34 


2002hw 


0.0163 


12.91 


10.01 


10.55 


10.67 


2002jy 


0.0201 


19.50 


10.11 


10.65 


10.78 


2003U 


0.0279 


17.93 


10.33 


10.48 


10.73 


2003cq 


0.0337 


25.75 


11.23 


11.34 


11.53 


2003ic 


0.0542 


54.79 


11.83 


11.97 


12.03 


2004L 


0.0334 


13.39 


10.50 


10.62 


10.80 


2005eq 


0.0284 


17.15 


10.82 


10.94 


11.12 


2005hc 


0.0450 


12.42 


10.45 


10.55 


10.67 


2005hj 


0.0570 


4.17 


9.50 


9.63 


9.77 


2005ir 


0.0753 


8.38 


10.16 


10.24 


10.33 


2005ki 


0.0208 


26.98 


10.86 


10.95 


11.50 


2005mc 


0.0261 


18.58 


10.87 


10.97 


11.26 


2005ms 


0.0259 


11.37 


10.24 


10.32 


10.65 


2006S 


0.0329 


15.05 


10.45 


10.60 


10.76 


2006ac 


0.0236 


16.43 


10.86 


11.11 


11.34 


2006ak 


0.0391 


10.34 


10.69 


10.78 


10.98 


2006al 


0.0691 


5.63 


10.27 


10.35 


10.45 


2006an 


0.0651 


2.62 


7.39 


7.68 


7.94 


2006ar 


0.0229 


6.73 


9.93 


10.22 


10.35 


2006az 


0.0316 


18.46 


11.48 


11.67 


11.75 


2006bd 


0.0265 


21.49 


10.95 


11.23 


11.32 


2006bq 


0.0215 


12.77 


11.00 


11.10 


11.38 


2006bt 


0.0325 


16.11 


11.23 


11.43 


11.52 


2006bz 


0.0277 


7.24 


10.27 


10.56 


10.65 


2006cc 


0.0328 


12.09 


10.83 


11.04 


11.12 


2006cf 


0.0423 


17.97 


10.92 


11.02 


11.18 


2006cg 


0.0288 


5.43 


9.94 


10.06 


10.33 


2006cj 


0.0684 


9.64 


10.57 


10.65 


10.75 


2006cp 


0.0233 


15.15 


10.07 


10.39 


10.45 


2006cq 


0.0491 


20.87 


10.90 


11.03 


11.17 


2006cs 


0.0241 


10.17 


10.89 


11.23 


11.28 


2006ef 


0.0170 


8.95 


10.24 


10.31 


10.87 


2006cj 


0.0192 


6.25 


10.83 


11.16 


11.22 


2006nz 


0.0372 


2.90 


10.03 


10.12 


10.31 


2006oa 


0.0589 


8.89 


8.80 


8.96 


9.09 


2006ob 


0.0582 


28.33 


11.14 


11.24 


11.35 


2006on 


0.0688 


4.95 


10.14 


10.23 


10.35 


2006or 


0.0217 


17.90 


11.01 


11.21 


11.43 


2006tc 


0.0321 


12.48 


10.43 


10.55 


10.74 


2007F 


0.0242 


12.22 


9.88 


9.97 


10.26 


2007R 


0.0296 


16.06 


11.17 


11.31 


11.47 


2007ap 


0.0162 


7.55 


10.24 


10.32 


10.87 


2007ar 


0.0534 


25.75 


11.61 


11.65 


11.69 


2007ba 


0.0391 


18.99 


11.09 


11.19 


11.37 


2007bc 


0.0219 


14.86 


10.83 


10.94 


11.24 


2007bd 


0.0320 


14.01 


10.68 


10.77 


11.01 


2007bz 


0.0227 


5.66 


9.45 


9.53 


9.84 


2007ci 


0.0191 


15.42 


10.63 


10.72 


11.28 


2008L 


0.0187 


7.32 


10.18 


10.27 


10.82 


2008af 


0.0346 


29.23 


11.53 


11.76 


11.81 



SN la Hosts and Calibration 

TABLE 12 — Continued 



SN 


Z C'MB 


90% i-band radius 


M*- M* M* + 






(kpc) 


(log M ) 


2008bf 


0.0253 


20.39 


11.25 11.59 11.64 



